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Abstract 

We consider the minimal supersymmetric extension of the Standard Model allowing both 
Dirac and Majorana gauginos. The Dirac masses are obtained by pairing up extra chiral 
multiplets: a singlet S for U{1)y, a triplet T for SU{2)^ and an octet Og for SU{3)c with the 
respective gauginos. The electroweak symmetry breaking sector is modified by the couplings 
of the new fields S and T to the Higgs doublets. We discuss two limits: i) both the adjoint 
scalars are decoupled with the main effect being the modification of the Higgs quartic coupling; 
ii) the singlet remaining light, and due to its direct coupling to sfermions, providing a new 
contribution to the soft masses and inducing new decay/production channels. We discuss the 
LSP in this scenario; after mentioning the possibility that it may be a Dirac gravitino, we 
focus on the case where it is identified with the lightest neutralino, and exhibit particular 
values of the parameter space where the relic density is in agreement with WMAP data. This 
is illustrated for different scenarios where the LSP is either a bino (in which case it can be a 
Dirac fermion) or bino-higgsino/wino mixtures. We also point out in each case the peculiarity 
of the model with respect to dark matter detection experiments. 
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1 Introduction 



A remarkable fact of nature is that the hght fundamental fermions appear in the smallest represen- 
tations, the singlet or fundamental, of the Standard Model SU{3)c x SU{2)m x U{1)y symmetry 
group. Could larger representations be present at higher energies and be discovered at the LHC? 
Such particles are in fact predicted by low energy supersymmetry as gauginos, superpartners of 
the gauge vector bosons. These gauginos can appear as fermions of Majorana (with only two 
degrees of freedom) or Dirac (with four degrees of freedom) type. In the Minimal Supersymmetric 
extension of the Standard Model (MSSM) the gauginos are Majorana. Obtaining Dirac masses 
would require pairing them with additional Dirac Gaugino adjoint (henceforth DG-adjoint) states: 
a singlet S for C/(l)y, a triplet T for 5C/(2)^ and an octet Og for 5f7(3)c. 

Construction of models with spontaneous breaking of supersymmetry may lead to non-minimal 
extensions of the Standard Model that include Dirac gauginos. The presence of the required 
DG-adjoints in the light spectrum is often motivated by the presence of an underlying N = 2 
supersymmetry, that pairs them with the vector multiplets. Such a scenario was first introduced 
by Fayet [1] as a way to give masses for gluinos while preserving R-symmetr}!^. The soft nature 
of these masses requires a modification of the interaction as was shown using D-term breaking 
by [21 [3]. More recently, Dirac gauginos have arisen in models with an extra dimension where 
supersymmetry is broken by a Scherk-Schwarz mechanism (see for example [H [5] ) . More precisely, 
they are a combination of two Majorana fermions with mass given by 1/2R (half of the compact- 
ification scale 1 /R) , one given by the (mass shifted) massless mode and the other from the (mass 
shifted) first Kaluza-Klein state, thus, the DG-adjoints originate as "half of the first Kaluza-Klein 
excitation" . It was noted there that the soft masses are UV finite and do not exhibit the usual log- 
arithmic sensitivity to the UV cut-off [6]. This property was shown to persist in four dimensional 
models as being peculiar to the Dirac nature of the gaugino masses, and denoted as supersoft 
in [71 El El [TOl [TT] where the important phenomenological implications of D-term supersymmetry 
breaking were first outlined. They were further studied in constructions of non-supersymmetric 
intersecting brane models [12\ [T3l I14j . More recent examples arise from the possibility of using 
calculable i?-symmetric models [31 HSl [161 El O US]. For instance, it was pointed out in [16] 
that they lessen the flavor problem in supersymmetric theories. This renewal of interest in such 
models has been motivated by the work of ISS [20]. Furthermore, deforming the ISS model with 
explicit breaking of ii-symmetry could leave states in adjoint representations [21j and allow the 
simultaneous presence of both Majorana and Dirac masses for the gauginos. The generation of 
Dirac gaugino masses can also be included [22] in the framework of "general gauge mediation" 

^ In [T] 7?-symmetry was later broken by Majorana masses for the DG-adjoints in order to avoid tachyonic masses 
for their scalar components. An issue that was recently solved in [22] and independently for an explicit model in 
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[231 [Ml |25l ESI [271128]. 

In this work, we will study a minimal extension of the MSSM that incorporates both Majorana 
and Dirac gaugino masses. The field content is that of the MSSM, supplemented with the DG- 
adjoints. The MSSM renormalisable Lagrangian is then supplemented by i) the DG-adjoint kinetic 
and mass terms, ii) the Dirac gaugino masses, iii) coupling of the singlet S and the triplet T to 
the Higgs doublets with strength and At respectively, iv) the DG-adjoint scalar soft masses 
and trilinear terms. 

For phenomenological issues the strength of the coupling of the Higgs to the DG-adjoint is of 
particular importance. In general such couplings are arbitrary and subject to diverse phenomeno- 
logical bounds as discussed in [8]. Inspired by extra dimensional models, and in order to make the 
role of = 2 manifest, one can assume that the two MSSM Higgs doublets originate from a single 
hypermultiplet of an underlying N = 2 supersymmetry. These models with combined N = 2 and 
N = 1 sectors were introduced in [HI [29]. In this case, the couplings As and At are related to 
the gauge couplings by = 2 supersymmetry. In this work, we will arbitrarily take the values of 
both couplings to go from zero to their tree level N = 2 value to illustrate the model. 

Some particular signatures of these models at collider experiments have been studied in [301 
[3T| [32[ l33l l34[ [35] . They stressed that the Dirac/Majorana nature of the gluinos affects the 
distribution of produced squark states. It was pointed out in [32] that the pair creation of scalar 
octets at the LHC will have a peculiar signature through cascade decays giving rise to a burst of 
eight or more jets together with four LSPs as well as through a resonance due to the decays into 
gluons or a tt pair at the one-loop level. 

In this work we will focus mainly on the fate of dark matter in this class of models. It is by 
now an important quality of the i?-parity preserving versions of the MSSM that they provide a 
natural candidate for Dark Matter, the lightest supersymmetric particle (LSP). A particular case 
is when the LSP is identified with the lightest neutralino. While in the MSSM the latter is a linear 
combination of the four neutral fermions, given by the bino, the wino and the two Higgsinos, it 
is now a linear combination of six states, the singlet and triplet fermions adding to the previous 
four. In this work we shall not try to give an exhaustive discussion of such a situation but try to 
answer such questions as: Are there parameter regions where the neutralino is a good dark matter 
candidate? How does the situation compare to the MSSM? An early study [36] of dark matter in 
a related model to that considered here focussed on the bino LSP, for a very particular case with 
vanishing Xs and At, and assumed dominance of the exchange of sfermions and thus neglecting, 
for example, the exchange of Higgs or gauge bosons. It concluded that the bino annihilation cross 
section can be enhanced which might help to obtain a smaller relic abundance than in the MSSM. 

The present paper is organized as follows: Section 2 presents the model and defines our 
notations and conventions. Section 3 gives a comprehensive discussion, which we believe to be 
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missing in the present literature, of the electroweak breaking sector. The discussion fohows the 
same hues as in |14] which specializes to a model with combined N = 2 and = 1 sectors in the 
limit of very large soft masses for the scalar components of T and S. It differs from the usual 
extension of the MSSM by the couplings to the Higgs of a singlet, whose vacuum expectation value 
(vev) is not necessarly related to the supersymmetric Higgs mass term ^u, and a triplet (see [3711381 
EHl HOI IHlllS] ) • While the triplet scalar is required to be heavy by electroweak precision tests, the 
singlet can be either very heavy and integrated out or remain light with sensible mixing with the 
ordinary MSSM Higgs states. We discuss both limits. In section 5, we briefly discuss the gravitino 
LSP and then focus on the case of a neutralino LSP. The corresponding mass matrix is exhibited 
and the nature of the lightest eiganstates is studied for some particular limits, in particular the 
necessary condition for having a Dirac fermion LSP is given. Section 6 presents numerical results 
for the relic adundance and the corresponding signature at direct /indirect detection experiments 
are stressed. 

2 The model 

The particle content of the model is presented in table [TJ The MSSM matter fields acquire masses 
through the Yukawa superpotential: 

WYukan,a = yJ/u^Qj • " y^^d?Qj • Hd - y>?Lj • Hd (2.1) 

and the usual soft breaking terms: 

+^^u'iQj ■ Hu - A'idtQ, ■ Hd - Ai^e^Lj ■ Hd + c.c. (2.2) 

where the bold characters denote superfields. Here, i,j are family indices and run from 1 to 3. 
The 3 X 3 y matrices are the Yukawa couplings. The "•" denotes SU (2) invariant couplings, for 
example: Q ■ Hu = ulH^ — diH:^. 

In order to have Dirac masses for the gauginos, additional fields in the adjoint representations, 
the "DG-adjoints" , are introduced. We define the superfields: 

s = s + V2exs + --- (2.3) 

T = T + V2exT + --- (2.4) 
Og = Og + V2exg + --- (2.5) 
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Names 




Spin 


Spin 1/2 


Spin 1 


SU{3), SU{2), U{1)y 


Quarks 
(x3 families) 


Q 


Q = {uL,dL) 


{uL,dL) 




(3, 2, 1/6) 
(3, 1, -2/3) 
(3, 1, 1/3) 


Leptons 
(x3 families) 


L 


{VeL,eL) 


iyeL-,(^L) 




(1, 2, -1/2) 
(1, 1, 1) 


Hiffgs 


Hu 
Hd 


\ U ' U/ 


\ u ' n/ 




(1, 2, 1/2) 
(1, 2, -1/2) 


Gluons 






[= 9a] 




(8, 1, 0) 


W 






r r -T ~ 1 TT-rOn 




(1, 3, 0) 


B 






N ^] 


B 


(1, 1, ) 


DG-octet 










(8, 1, 0) 


DG-tripIet 


T 


{ro,r±} 
N {So^,s±}] 


[= {W'^,W'^}] 




(1,3, ) 


DG-singlet 


S 


s 


XS 




(1, 1, ) 



Table 1: Chiral and gauge multiplet fields in the model 



where S = -j^{Sr + iSj) is a singlet and T = Ylia=i 2 3 ^^"^ SU{2) triplet parametrized as: 



i 2 ' ^ 2 ' "2 

1/ To ^T+\ 
2 VV2T_ -To y ' 

ro = 4s(rR + ir,), T+ = ^{T+n + iT+i), = ^{T^R + iT.j), (2.6) 



and fj'^ are the Pauli matrices. Their quantum numbers are presented in Table [H 
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Due to the presence of these extra fields, the gauge kinetic terms are modified to become: 

>C„ = J d^xSe [ ^Mi Wf Wi, + ^M2tr(W?W2a) + ^M3tr(W^W3a) 
+V2m?oWiaS + 2V2m^^tr(W2aT) + 2V2m^^tr(W3aOg) ] 
+ fct^xcfecfe (^*|e*^-i*i + /i.c.) (2.7) 

where Vj are the vector and Wj^, the corresponding field strength supcrfields associated to U{1)y, 
SU (2) and SU (3) for j = 1,2,3 respectively. Here, we have introduced spurion superfields to take 
into account the generation of gaugino masses: 

Mi = l + 2eeMi (2.8) 
maiD = OaVniD (2.9) 

The Dirac gaugino spurion superfield can be written as = —^DDDa^i- This mass originates 
as a D-term if Xj is identified as a vector field Xj = V'/Mj, or as non vanishing F-term by writing 
Xj = 2S^5]/Mj with S = OOF, where Mj is the appropriate super symmetry breaking mediation 
scale. 

The DG-adjoints may also modify the Higgs superpotential, since new relevant and marginal 
operators are now allowed: 



j d'^x(fO /iHu • Hd + + AsSHd • Hu + Mrtr(TT) + 2ArHd • THu (2.10) 

with the definition i7„ • iJ^ = H+H' - H^H^. 

Finally, the soft supersymmetry breaking terms for the scalars are: 

-/^Csoft = m]jJ\Hu\^ +m]j^\Hd\^ + B^{Hu- Ha + h.c.) 

+m%\S\^ + \bs{S'^ + h.c.) + 2m^tr(rtr) + BT{tr{TT) + /i.e.) 

+As\siSHd ■ Hu + h.c.) + 2AT\T{Hd ■ THu + h.c.) (2.11) 

Note that we did not include in the superpotential a cubic term tr(TTT) as this identically 
vanishes. Neither did we include linear and cubic terms in the singlet. The latter is due to the 
fact that we assume that the DG-adjoint appears due to some underlying N = 2 supersymmetry 
that forbids these terms. Of course a microscopic model explaining the origin of the soft terms 
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should also address the fact that the super symmetric Mg and Mt are assumed to take values of 
order of the electroweak scale, introducing an issue of scale hierarchy as does the Higgs /i-term. 

Below, we will give special attention to the scenario where the DG-states arise as a result of an 
N = 2 extension of the gauge sector. In this case, if the Higgs multiplets and are assumed 
to form an = 2 hypermultiplet then A5 and \t are related to the gauge couplings, at the N = 2 
scale, by: 

\s = V2g']^, At = V2g^, (2.12) 

where g' and g are the U{1)y and SU{2) gauge couplings respectively. The factor 1/2 in Xg arises 
from the U{1)y charge of the Higgs doublets. 

3 Electroweak scalar potential 

We turn now to the electroweak scalar potential. This receives contributions from three sources: 

VeW = y gauge + Vw + ^so/t (3.1) 

The first is a contribution from the gauge kinetic term (j2.7p . Integration on the spinor coordinates, 
and going on-shell, leads to the U{\)y and SU{2) D-terms: 

Di = -2miDSR + DP with Df> = -g' ^Yj^*ipj (3.2) 

j 

= -^/2m2D(T" + r'^t)+ 2)^(0) ^ith D^(0) = _5^c^*^(^^. (3.3) 

j 

where ^pj are the scalar components of matter chiral superfields, whereas ^ and D'^^^ are the 
D-terms in absence of Dirac masses. The resulting Lagrangian contains terms of the form: 

C-gauge ^ -miD^XSa " m2Z)tr(Af XTa) " \dI - ^tr {D^D^) (3.4) 

where we can identify the Dirac components of the gauginos Ai = B' and = W"^ as given in 
Tabled! 

The contribution from the DG-triplet to is: 

1 M|r_|2-|r+|2) ^{T+T*-T,Ti)\ 
' 2 \V2{ToTi - T.T^) -(|r_|2 - |r+|2) ) 
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which vanishes in electrically neutral vacuum, where: 

< T+ > = < T_ >=< H+ >=< i7_ >= 0. (3.6) 
The contribution to the scalar potential of the neutral fields is then given by: 

Vgauge = ^mj^Sl - 2miDSRDP + ]^dP^ + 2mi^r| - 2m2D7ijZ)f + ^Z)f' (3.7) 

where we have dropped the generator label, D2 = ^f'^ ~ ^' only non-vanishing compo- 

nent. 

The second contribution comes from the superpotential (j2.10p : 

W = (-A^ + A5S)(H°H°~H+H^) + ^S2 + ^(T°T° + 2T+T-) 

-At(H^T°H° + H^T°H;) - ^/2At(H-T+H° - H^Jt-H^). (3.8) 
Keeping only the neutral components, it reads: 

Vw = \MsS + XsH^'^H^S + IMtT" - At^XI' + " >^sS + \tT'^\\\H^^'' + l^op) (3.9) 

The third source Vsoft is due to soft supersymmetry breaking terms for the scalars ()2.1ip . We 
define: 

Hu- ^ , H,- — (3.10) 

then, all together, the scalar Lagrangian for the neutral fields is now given by: 

Vew = {rnl^+^^')\H'u\' + {rnl,+Mf-B,{^^^^^^ 

+(A| + A|)|i/°i/°|2 (3.11) 
+^(M| + m| + 4m?^ + Bs)Sl + ^(M| + m| - Bs)S] 

+^(M| + + Aml^ + Bt)T^ + ^(M| + m| - BT)Tf 



+ 



■^{Sl + Sj) + -^{Tf + T|) - V2f,iXsSR - XtTr) - XsXt{SiTj + SrTr) 
X ri ij0|2 , I H-0121 



Note that the MSSM potential is given by the first line. All the parameters are chosen to be 
real. We are left with four neutral fields Sr, Si^Tr^Tj satisfying equations of type: 



a = R,I for {x = S,y = T) or {x = T,y = R) 



(3.12) 



with solutions of the form: 

_ V^gM^a - VygXsT 



xa ya 



a = R,I for [x = S,y = T) or (x = T,y = R) 



(3.13) 



where: 



while 



and: 



mIj 



Ml + ml + Amjn + Bs + A|(|i7°p + \H°af) 

Ml + ml-Bs + \UK\' + \H'd\') 



t0\2\ 



XsT = -XsXt{\HI\'' + \H^^ 



t0\2\ 



(3.14) 
(3.15) 
(3.16) 
(3.17) 

(3.18) 



VsR = V2^lXs{\H°J' + \H'af)-9'mlD{\H^f-\H'df) 
-^{Ms + As){H',rH'^^ - H^,jHl) 

Vtr = -V2^,XT{\H^J' + \H',\') + gmMW - 

Vsi = -^{Ms-As){H^aRHli + H^diHlR) 
Vti = +^{Mt - At){H^,j,hIj + H^,jHl^) 



(3.19) 

(3.20) 
(3.21) 
(3.22) 



We are not going to pursue exact computations. Instead, we will consider the case with 
m|, m\ S> m\ in which case the formulae simplify as we can neglect all terms proportional to 
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+ in particular Xst- This gives Xa ^ Vxa/M^^, i.e. 

As, 



Ml + m% + 4m2„ + Bs 



Tr ~ 



Si ^ -\s 



{Ms-As){H^,rHI, + H^,,HI^) 
V2(M| + m| + Aml^ + Bs) 



We are interested by the case of CP neutral vacuum, i.e. H^^ = H^j = which implies 
Sj = Tj = 0. From now on, we will drop the indices and defin^: 

< H\ >=< hu> = Vu = vsf3, < H^R >=< hd >= Vd = vcp, 0^/3^1 

<Sr> = Vs <TR>=vt (3.24) 

where we denote: 

cp = cos /5, S/3 = sin/3, tp = ian(3 

C2/3 = cos 2/3, S2/3 = sin2/3 (3.25) 



vt 



2(MJ + m| + 4mf^ + Bs) 



2(M2 + m|, + Amlj^ + Bt) 



g miDC2i3 + V2nXs - -^{Ms + As)s2p 
-gm2DC2f3 - V2^i\t + "^(^^ + ^t)s2/3 



. (3.26) 



Electroweak precision data give strong bounds on the expectation value of the DG-triplet as it 
contributes to /) ~ 1 + aT = 1.00041°;°°°^ |i3]. Thus we require: 



A/9 =i 4^ < 8 • 10- 



(3.27) 



which is satisfied for vt ^ 2> GeV. Here we will allow the Dirac and Majorana masses to vary 
arbitrarily and will satisfy this bound by taking rriT large enough. For instance, for Ms,As,^ ~ 
200 GeV and all couplings of order one, this is satisfied for niT ^ 1 TeV. 



^With this convention v ~ 246 GeV, 'g')^+g% ^ = M|. 
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Integrating out the DG-adjoints S and T, leads then to the effective tree-level scalar potential 
at first order in Xs,t 



Vew 



32 



\2 I \2 

1 



1 - hi) - V2^i\s{hl + hi) + ^/2Ag(Mg + 

1 [-ffm2g(^^ - hi) + V2fiXT{hl + /i^) - V2Xt{Mt + At)/^^^^]' 
8 M|, + m2, + 4m2r, + Sr 



(3.28) 



This decomposes into three parts: 



Vew = Vo + Vi + V2 



The first part: 



Vn 



h, + 



2 " 2 

is the MSSM contribution. The second, 

Vi = 



'hi - B^huhd + 



2\2 



32 



iK - K) 



X2 _L X2 

^5 "I" t.2 ^2 

^ "■u'^d 



(3.29) 



(3.30) 



(3.31) 



is a quartic term. The third contains the explicit dependence on the mass parameters of the 
DG-adjoints. We will illustrate this in taking a few limits, keeping the other parameters fixed: 



One limit is to take Mg — oo and Mt — do. In this case: 



V2 



'^S + ■^T 1.2 1.2 



(3.32) 



meaning that the complete DG-adjoint supermultiplets have been decoupled and one is left 
with the MSSM electroweak scalar potential. 

A second limit is to take A5 — >^ and At — switching off the superpotential couplings 
Vi — ^ (3.33) 



V2 
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n 



00 V , 



32 4mfr, -^m^^Ml^Bs 32 Mm^^ + ml + + Bt 



2\2 



K - hi) 
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This shows that the effect of Dirac masses is to decrease the quartic couphng, and make it 
vanish in the absence of other masses. 

• A third one is to take niiD oo and m2D oo 

V2 -'^^{hl-hlf (3.34) 

the Higgs quartic term of the MSSM due to D-term is cancehed and all the quartic couplings 
are generated by superpotential coupling between the Higgses and the DG-adjoints. 

• As a fourth one, we consider the case of interest in the rest of this work: ms — > oo and 

rriT — > oo 

\2 _L \2 

V2 — > (3.35) 

showing that the MSSM scalar potential is supplemented with a quartic term that lifts the 
D-term flat direction = H^- 

3.1 Integrating out the adjoints 

In this last limit, the vacuum expectation values of the DG-adjoints can be neglected, allowing to 
write an SU (2) invariant effective scalar potential. It can be put in the usual form |44l H5] which 
parametrizes the two-doublet potential: 



+1 [\{9' + 9") + Ai] + I [\i9' + 9") + A2] i\H^\'f 

+ [\i9' - 5") + A3] l^rfPl^.l' + [ - + (Hd ■ H^){Hl ■ HI) 

+ {^{Hd ■ Huf + [Ael^dp + Xr\H^f] {Ha • F„) + h.c.) (3.36) 



where now: 



A3 — 2A2^ A4 — A^ — 

Ai = A2 = A5 = A6 = A7 = 0. (3.37) 
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As in the MSSM we expect sizeable one-loop corrections to the above potential, but we will also 
have new and particularly important corrections due to taking m^, rriT large. We can approximate 
these new contributions by taking the leading logarithmic contributions from scalar loops involving 
only quartic vertices; the remaining diagrams involving scalars and sfermions will have the effect 
of removing the dependence on the cutoff. We can reasonably approximate this behaviour by 
replacing the renormalisation scale in the logarithms with v (a more precise calculation involves 
the full fermion mass matrices). We should then also evaluate the coupling constants at this scale. 
We can then write 

and can neglect the contributions to the other factors as subleading. For simplicity we have 
assumed Bs, Bt ^ rog, m^; this is not valid in the scenario of for example [22] where they are of 
equal order in magnitude, both being generated at one loop. 

Note that there is also a contribution to the Higgs mass parameters proportional to 

A| J, log ^ ^ , which is absorbed into the renormalisation of m?jj^,m?^_^. If these are large 

then we have a large fine tuning. 

3.2 Singlet Extension to the Higgs Sector 

As can be seen from equations ()3.2|3.3p . the sfermions obtain a mass proportional to the expecta- 
tion values of Sr^Tr. Of these, we can ignore the Tr contribution due to the strong constraints 
upon it as discussed above, but {Sn) may be non-negligible. Together with the dominant one loop 
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effect (containing gauginos or Sn,Tji scalars) we find 

ml = (n^f + 2m,n{Sn)9% + Y^ig' f^ log ^ + log ^ (3.39) 



where, since the action that we are using is an effective one, we have included the soft masses 
induced by tlic supersymmetry breaking sector (via gauge mediation or otlicrwise) otlier tlian 

through the Dirac gauginos as rn^^^ For example, in gauge mediation these are generated at 
two loops; the gaugino masses miD,m2D,rnx^,mx2 (the latter two are the total gaugino mass, 
including both Dirac and any Majorana effects, defined as the location of the pole in the (AA) 
propagator) are generated at one. The one loop fluctuations around our effective action then 
actually appear in gauge mediation at three loops. 

The term proportional to (Sr) clearly gives a negative contribution to negatively charged 
states, so we must ensure that this does not dominate. This is particularly important for models 

where m^'^^ vanishes, which can occur for certain gauge mediation models with purely Dirac 
gauginos (in such cases MsjMt = 0). One way to ensure positivity is to take a large ms,mT, 
which reduces {Sr) and increases the loop effects. We then have (where j,^ = m| j, + Bs,t) 



This is easy to arrange for large msf^ , for example for small A5 or large miDC2f} <C ms we require 
(from SU (2) singlets, which provide the strictest constraint) 

and therefore for m\j^ ~ we have msj^ ^ 3f + 0.94(mAi — v) + ... > 750GeV, and for mfjj = 



xm|^ we require m^^ ^ t^v/ a/— (2x + l72)logx. If we consider for example that the masses are 
generated in gauge mediation, then rriiD ~ Xxg'^, ms ~ Bs ~ ^x^^ ^'^^ ^ messenger coupling 
Ax and effective supersymmetry breaking scale A; this gives x = (5')^ ~ 1/8.4 and thus mg^ > 
620GeV. However, this limit is only valid for A5 ^ g'/V2, as the gaugino masses are smaller than 
V. For A5 ~ g' /V2 we can place a bound by setting ^ ^ v: 

(5 ) A 9 log a; > (3.41) 



2m5^(4a; + l) ' 47r2 
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1/3 



and thus ms„>v 



27rV(\/x(l + 4x)logx) 



> 650Gey. 



Now consider when ms is smah, Xs{Sr) = ft is large and adds to or replaces the n term 
while AsXs{Sfi) contributes to S/i. The loop terms are no longer significant, and we must simply 
impose that 



This will be in general quite a stringent constraint on the parameter space of the microscopic 
model. 

Keeping S in the light spectrum results in the effective potential (j3.12p with T = 0. The 
MSSM spectrum is extended by one CP even Higgs and one CP odd state, corresponding to Sr 
and Si respectively. This falls in the class of models studied by [46j . In addition to the possibility 
discussed there of new decays of the lightest Higgs into pairs of CP odd states (as can be seen from 
Eq. I3.12P we would like to stress the new feature that the Higgs has different decay/production 
channels due to the mixing with Sr which couples to the D term as 



The effect of this at colliders is very model dependent. It is important, for light sfermions, when 
msji is smaller than or comparable to rrin). For example, if sneutrinos are arranged to be light 
then the Higgs may decay to them and then to neutrinos, plus the LSP. 

4 Fermion mass matrix 
4.1 The neutralinos 

There are six neutral fermions of interest: the higgsinos and H^, the gauginos, bino B and 

wino and the DG-adjoint fermions B' and wino W'^. The mass terms for these fields have 
different origins: 

• bino, wino and DG-adjoints Majorana masses: 



{m\ff + 2miD-^g'Yi > 0. 



(3.42) 




(3.43) 




(4.1) 



• bino and winos Dirac masses: 



m2DW"W' 



-la 



m\£)BB' + h.c. 



(4.2) 
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• The gauge interactions between the gauginos, the higgsinos and the scalar Higgs: 



(4.3) 



leading to 



mz 



sin 0w{s,3 HlB - iilB) + cos dw[c^ &dW^ - S/3 + h.c. (4.4) 



• The superpotential (|2.1U|) leads to couplings between the DG-adjoint fermions, Higgs and 
Higgsinos 



\s [Rd ■ HuB' - Hu ■ HdB') - At ■ {a'Hd)W" + Hd ■ {a'Hu)W 



(4.5) 



giving 



mz 



9' 



(4.6) 



• The fi term in the superpotential W contributes to the higgsinos masses 



(4.7) 



All the previous terms together describe the resulting mass matrices for both neutral gauginos 
and higgsinos when both Majorana and Dirac term are present. The neutralino mass matrix A^Oj 

in the {&, B, W"^, V^o, H^, Hi) basis is: 



/ M[ miD 

miD Ml 



m2D 

^^Y^mzswSfs -mzswcp -^^^rnzcwsp mzcwcp 

\ ^^^mzswcp mzswsp 



m2D 

M2 



^^Y^mzswsp ^^^^mzswcp 

-mzswCf3 mzswSjs 

-^^^mzcwsp -^^^mzcwCfs 

mzcwcp -mzcwsis 

-12 



\ 



^^^mzcwcp -mzcwS/3 



-1^ 







(4.8) 

where cw = cosOw, sw = sin6'vy. This may be diagonalised by a unitary matrix A'' such that 
^diag ^ N*MoN^; perturbative expansions for N in various limits are given in appendices Rl 
andE 
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4.2 The charginos 



The chargino mass matrix describes the mixing between the charged higgsinos , and the 

charged gauginos W"^ and . The corresponding mass terms have the same origin as those 
presented in section WA\ There are Dirac masses: 

- - m2DW+W'- + h.c. (4.9) 

There is the usual ii-teim 

- I,H+ ■ H, + h.c. (4.10) 
and finahy there are mixing terms of the gauginos with the Higgsinos: 

- V2mw sin PH:^W~ - V2mw cos + h.c. (4.11) 
and with the DG-adjoints 

2A 2A 

- -^mwcpH^W'- + -^mwspHjW'^ + h.c. (4.12) 

The mass terms for the charginos can be expressed in the form 

- ^{{v-fMchV^ + {v+fMl,,^v- + h.c) (4.13) 

where we have adopted the basis = (W'~^ ,W~^ , H:^), v~ = {W'~ ,W~ ,H^). Collecting all the 
terms presented in equations (|4.9p - (l4.10p . (14. 2p and (I4.12P leads to the chargino mass matrix : 



M. 



Ch 



^ M'2 m2D "^mwcp 

m2D M2 V2mwsi3 \ ■ (4.14) 



This nonsymmetric matrix is diagonalized by separate unitary transformations in the basis 
and M'^i^^ = WMchV, where U and V are unitary. 

5 The LSP dark matter 

The model has i?-parity so the LSP is stable. Here, we assume the LSP to be the (lightest) 
neutralino and we would like to study its relic density and see how it fits with actual bounds 
from WMAP. However, before doing so we would like to make some comments on the fate of 
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gravitinos in this scenario. Gravitinos play two major roles in our model: 1) there is the issue of 
i?-symmetry breaking by the gravitino mass m-^j^ and inducing Major ana masses for the gauginos 
2) the gravitino may be the LSP for instance in models of gauge mediation. Depending on its 
mass it can also change the LSP relic abundance. 

First, let us discuss the first issue. In AT = 1 supergravity the gravitino mass needs to be 
non-zero to cancel the cosmological constant after supersymmetry breaking. It is proportional to 
the vev of the superpotential and thus breaks i?-symmetry. If one insists on avoiding the breaking 
of R symmetry one option is to enhance the gravitational sector toiV = 2(=Ar = li©iV=l2) 
supersymmetry. The supersymmetry breaking preserves i2-symmetry when the two gravitino 
masses are equal. The set up is then a Dirac gravitino made of two gravitinos '^I'^j^ i and -03/2 2 
with the same mass m3/2, each of them coupling to a different sector with N = \\ for the first 
and iV = I2 for the second. The coupling strength is in both cases given by but may 

couple different gravitinos to different sectors. We will not present here an explicit realisation as 
it goes beyond the scope of this paper. 

The second issue is very model dependent. A gravitino LSP can be produced in two different 
ways: either through thermal production or through decays of unstable sparticles. The relic 
density is expected to receive contributions from the two processes by an amount that depends on 
the particular spectrum and in the specific thermal history of the universe (for instance depending 
on the value of the reheat temperature), the possible decays of the inflatons, moduli, etc. The 
addition of a second gravitino might help the gravitino LSP to be a dark matter candidate by 
increasing the number of relativistic degrees of freedom at the freezout temperature. However, 
with a second sector, it will make the story even more dependent of the details. We will not 
pursue this scenario further here. 

The neutralino sector that we shall consider differs from that of the MSSM in two aspects: 

• The LSP is now a linear combination of six states. There are now two additional states 
compared to the MSSM case. 

• In addition to the MSSM parameters, we have six additional parameters: 

M{,M2,mi£),m2_D, As and At (5.1) 

In this work, our purpose is to find a region of these parameters where the relic density is 
compatible with the dark matter made of a thermally produced neutralino LSP and compare this 
with the case of MSSM. We will assume both the S and T scalar are very heavy and decoupled 
from the thermal bath. 
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5.1 The Dirac case 



The first simplest case we will consider is the LSP to be mainly bino like. 



/ 




/2A 

\ ^mzswcp 



miD 




-mzswcp 
mzswsis 







m2D 
-^^^Y^mzcwsp 
-^^^mzcwCfs 



With /X > miD The DG-adjoint vevs are now 






m2D 


mzcwcp 
-mzcwsis 



^^Y^mzswsp 

-mzswCf3 
-^^^mzcwsp 
mzcwCjs 




vt 



2{M§ + m| + 4mf^ + Bs) 



2(M'2 + m^ + Aml^ + Bt) 



g miDC2i3 



"^mzswcfs \ 

mzswSjB 
-^^^mzcwCfs 
-mzcwSfB 

J 

(5.2) 



gm2DC2i3 - \/2/iAr + ^i^T + ^t)s2/3 



(5.3) 



To have a pure Dirac LSP, the lightest two eigenvalues of the neutralino mass matrix [5^2] most 
form a pair of equal magnitude but opposite sign. It is straightforward to show that in the case 
that both Xs,Xt take their N = 2 values we will only have Dirac neutralinos. To do this, we 
note that solving for the eigenvalues of the neutralino mass matrix Ai means solving the equation 

/(A) = det(A^ — A) = 0. For purely Dirac states, /(A) = n?=i(-^^ ~ ) eigenvalues ±aj, which 
is true if and only if /(—A) = /(A). By examining the coefficients of A^ and A^ ( that of A^ being 
automatically zero due to the tracelessness of A^) we find that this requires 



A^ : 2^cpspMl 



: 2ncpspMz 



2 2 / 1 2At\ 2 2/^1 2A5 



9 



2\i 



c^(i-^) + 4/1 1 



9^ 



2Xs_\ 



{9'? 
= 







(5.4) 



proving the above assertion. 

If we assume a mostly bino/f/(l) adjoint LSP, by assuming that miD < rn2D ^ ^ we can 
expand the LSP eigenvalues to next to leading order in m\^/ ii,m\j^/ ^ to find 



rriLSP = + 



Mis 



2 „2 



i9'y^ 



V2\sg'{sl - 4) ± (2A| - {g' f)cps/3 



+ 



(5.5) 
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Assuming that (2A|. — {g')'^)cpsp is negative, the LSP takes the upper sign, and the mass sphtting 
is given to lowest order by 



AniLSP 



,M|.^ (2A|-(gT) 
-2 r^TT^ CpSp. 



id') 



(5.6) 



Note that this reduces to the result of [36] when A5 — > 0. The eigenstates at this order can be 
read off from the rotation matrices, given in appendix El 

For the case of a Higgsino LSP, where rriiD ~ ^ we find 



rriLSP + 



/3 



4)M| 



gg'miDm2D 



Cwg'\TiniD + SyygXsm2D 



+ 



M| \'^cp^pMl{cl^g' \TrniD + g\srn2Dsly 



2 \2-i 



g'^{g')^mfjjm? 



At 



1D"''2D 



M',[ciy{gr{g' + 2A^)mf^ + + 2A^)m^^s^^ 



2c/2(g,')2m2 2 



=F 



g^{g')^mlj^mljj 



'■1D"''2D 

{g'fmlocW - 2X^T)+9^rn^2Dsly{{g' f 



2Xi 



(5.7) 



where we have had to expand to the second order to obtain a mass splitting. 

If we suppose that the adjoint couplings take the N = 2 values at some scale Mn=2 and we 
run A5,At down to the supersymmetry breaking scale M]\f=i (equal to ms or itit) then we can 
generate a mass splitting since the adjoints do not couple to the matter multiplets, and the Higgs' 
wavefunction renormalisation also contributes. To leading order we have 



2A| - {g'f 
2A^ - 



a/a 



167r2 

2g^ 
167r2 



3|yt|' + 3|yfe|2 + |y,|2-l0(5 



i\2 



log 



.Mn= 
■Mn- 



3|yt|' + 3|yb|2 + |2/,|2-4g2 



1 ,Mn=2. 
log! I 



(5.8) 



If we assume that yt ~ 1 and the other couplings much smaller, then we find for a Bino LSP that 
the mass difference should be 



ArriLSP ~ 2C/3S/3 



log(— ) 



O.lSGeV 



,TeV tf3 



(5.9) 



19 



where, obviously, smaller splitting can be seen to correspond to larger values of tan /5 and smaller 
ratio Mn=2/ M]\f=i. For instance, taking fi = ITeV, Mjv=2 ~ lO-'^^GeV, and tp = 50 leads to 



^misp ~ 4.4 




GeV ^ 90MeV 



(5.10) 



while for n = ITeV, Mjv=2 ~ lO'^GeV, and tp = 50, we obtain 



ArriLSP ~ 0.3 




GeV w 6MeV. 



(5.11) 



The latter case could be realized in low cut-off scale models (such as for large extra-dimensions). 
5.2 Dark matter relic abundance 

In the DG model, as in the MSSM, there are several scenarios that lead to a relic abundance of 
the neutralino LSP, Qh'^ ~ 0.11. These include naturally the MSSM- like scenarios 

• a mixed bino/Higgsino LSP that annihilates mainly into W pairs (or top pairs) 

• a mixed bino/wino or bino/wino/Higgsino that annnihilates mainly into W pairs 

• a bino that annihilates into fermion pairs when sleptons are light, a significant region of 
parameter space where this process occurs has been ruled out by LEP. 

• a bino that coannihilates with sfermions. 

• an almost pure bino with mass 2m^ ~ mh,A with efficient annihilation through Higgs 



New dark matter scenarios occur as well. Firstly in the special case of a pure Dirac bino LSP 
(or an almost pure Dirac bino), annihilation into light fermion pairs becomes efficient. This is 
because the process does not have a strong p-wave suppression as in the Majorana neutralino 
case [36j. Secondly the MSSM-like scenarios where B,W are replaced with B' and W can also 
occur. These are found either with Dirac masses or with Majorana masses. Typically, as we will 
see in case studies below, these involve more coannihilation processes. 

In the DG model, the dark matter detection properties can be quite different than in the usual 
MSSM. For one, the annihilation of a Dirac neutralino into light fermion pairs is not suppressed, 
even at u/c ~ 0.001. More generally, the content of the LSP that determines the coupling of the 
LSP to other particles can be different in the DG model. This is largely due to the additional 
B' or W components. For example, the spin independent elastic scattering rate on nucleons that 



exchange. 
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is dominated by Higgs exchange unless squarks are light, depends on the higgsino content of the 
LSP. In the DG model, the higgsino fraction is often suppressed leading to small rates. 

In the following subsections we present results for a few case studies. For each scenario, we 
find the parameter space that predicts a relic abundance compatible with the value determined by 
cosmological measurements, fi/i^ = 0.113 it 0.0034 [47j and examine the predictions for the direct 
and indirect detection rates. In all cases special attention is paid to the typical DM scenarios with 

either a B{B') or mixed B{B')/h or B{B')/W{W') scenarios (here h stands for H'^ or H^). The 
possibility of efficient annihilation through a Higgs resonance is a generic feature of all models 
where s-channel resonance can occur, we will however not consider this mechanism in detail as it 
requires fine-tuning of the model parameters. 

6 Results 

All numerical results are based on micrOMEGAs2.3 for the computation of the spectrum, the relic 
abundance, the elastic scattering rate as well as the annihilation rate o"f |o relevant for indirect 
detection [48\ I49j . We have implemented in this code the DG model described in section 3. 
The one- loop quark/squark corrections to the Higgs masses are computed as well as the A5, At 
corrections to the effective potential (I3.36p . The latter can increase the Higgs mass by a few 
GeV as compared with the MSSM. We have imposed the LEP bounds on charged sparticles as 
well as on the Higgs mass {mh > 111 GeV), allowing a large theoretical uncertainty for m/j since 
two-loops corrections are ignored. In all cases we fix tan (3 = 10, Af = — l.STeV and Mg = ITeV. 
Under these conditions the Higgs mass limit is easily satisfied. The elastic scattering rates are 
computed taking the micrOMEGAs default values for the quark coefficients in the nucleon [49j . 
Varying these coefficients could induce large corrections to the predicted rate. 

6.1 Pure Dirac masses : Mi = M[ = M2 = = 

In the case of pure Dirac masses Mid, M2D 0) the LSP can be a Dirac fermion provided Xs and 
Xt take their N=2 value, eq. ()2.12p . In general one expects a mass splitting generated by the N=2 
breaking effect, eq. (15. Qh . nevertheless it is possible to tune the parameters of the model such that 
the two lightest Majorana states are degenerate and make a Dirac fermion. The implications for 
dark matter detection of a Dirac fermion are important so this case is worth consideration. 
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6.1.1 Dirac fermion 



A Dirac neutralino, contrary to a Majorana neutralino, can annihilate into light fermion pairs with 
a large rate, thus offering an explanation to the excess of positrons seen by Pamela [501151] without 
spoiling the antiproton observations |52j [§. In particular the Dirac bino can have a much larger 
annihilation rate into leptons than into quarks when the mass of the right-handed sleptons are of 
the order of the bino mass. However, the Dirac neutralino has an effective vectorial interaction with 
quarks in the nucleon, this leads to potentially large rates for direct detection. For a Dirac fermion 
the spin-independent elastic scattering cross section receives dominant contributions from squark 
and Z exchange. To avoid exceeding the experimental bound \54\ [55] it is enough to 1) fix the 
mass of the squark that couples most strongly to the bino, the one with the largest hypercharge, to 
mji^ > 1. — 1.2 TeV and 2) suppress the higgsino LSP component such that the coupling of the LSP 
to the Z is reduced. To illustrate this we have computed the neutralino nucleon elastic scattering 
cross section in a scenario where M2D = l-5Mi/), /U = 1 TeV, Mj^ TeV and Mq^ = lorl.2 TeV. 

The value of the common mass for the right-handed sleptons, M^^, is adjusted for each value of 

Mid such that Qh'^ = 0.11. The lower value of the squark mass leads to a cross-section exceeding 
the CDMS limit for light neutralinos, see Fig. [TJ Note that one characteristic of scenarios with 
a Dirac particle as DM is that the spin independent amplitude for elastic scattering on nucleons 
can be different for protons and neutrons, for example cr^Op « c^o^ if Z exchange dominates. 

However the experimental limits on c-j^Op are extracted the amplitudes for protons and neutrons 
to be equal. In order to be able to compare directly with the experimental limit from CDMS |54j . 
we rescale the nucleon cross section and define an effective = {Zfp + (A — Z) fn)"^ / where 
fp{n) are the amplitudes on nucleons and Z = 32, A = 76 for Germanium. In this scenario, the 
rate for indirect detection is large, av\ii PS 2.5 X 10-26 cm^/sec for Mid = 300 GeV and the lepton 
channels are almost two orders of magnitude larger than the quark channels. 

For the remainder of this section we will consider the case where the LSP is a Majorana 
fermion, and the LSP-NLSP mass splitting in the MeV to GeV range as illustrated by eq. (j5.10p 
and eq. ()5.1ip . This results from shifts of the value of \s (say by only 1%) from the N=2 value 
in eq. ()2.12|) . Because of the small mass splitting between the LSP and the NLSP the X1X2 chan- 
nels will not suffer from a Boltzman suppression and will contribute significantly to the effective 
annihilation cross section that enters the standard computation of the relic abundance [48j . 

^Note that the recent Fermi results ^53[ on the total electron and positron spectrum is in agreement with the 
positron spectrum measured by PAMELA. 
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Figure 1: a) Effective neutralino-nucleon elastic scattering cross section as a function of m^o for 
a Dirac LSP case witli rriq = 1 TeV (blue) and rriq = 1.2 TeV(green). For each LSP mass, the 
common right-handed slepton mass is adjusted so that = 0.11. The CDMS limit is also 
displayed (black). 

6.1.2 Bino-LSP 

The main mechanism for annihilation of a bino LSP is through exchange of sfermions, the sfermions 
with largest hypercharge, the right-handed (RH) sleptons, giving the dominant contribution. In 
general this process gives only il/i^ ~ C'(l). This is because a oc m'io/m^- . and both the 

neutralino and the sfermion need to be near lOOGeV, that is near the LEP exclusion region, to 
reach fi/i^ = 0.1. Slepton coannihilation provides an alternative for reducing the relic abundance. 
With Dirac mass terms and nearly degenerate B and B', the processes X1X2 ~^ ff provide the 
dominant annihilation mechanism [36]. One can obtain Qh^ w 0.1 even with sleptons twice as 
heavy as the bino-LSP. This is displayed in Fig. [2] where we compare the slepton-LSP mass splitting 
that produces a relic abundance in the WMAP range in the DG model and in the MSSM. Here we 
show contours including a large theoretical uncertainty in the determination of Qh'^ = 0.11 it 0.026. 
In both cases we fix /i = 1 TeV and take M2D = 2Afiz)(DG model) or M2 = 2Mi(MSSM), so 
that m^o ~ Mid{Mi). For each neutralino mass we vary the RH slepton masses to find the given 
relic abundance contour. Because of the mixing in the stau sector, the stau turns out to be the 
lightest slepton. In this scenario the elastic scattering cross-section is small {c^p ~ 10^^^'^ pb) 
although within the reach of detectors such as Xenon [56j. The small cross sections can be linked 
to the LSP higgsino composition. The indirect detection cross section is also small in both models 
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Figure 2: a)Contours of Qh'^ = 0.094 (upper) and Qh'^ = 0.094 (lower) in the A(mf — rn^o) vs 
m^o plane in the MSSM(full blue line), the DG model(red dashed line) and the DG model with 
-^5 = (green dotted hue), b) Elastic scattering function of m^o - same colour code as 

a) 

ctvIq ~ 10~^^cm^sec. 

The prediction for the elastic scattering cross-section can be shifted significantly for a different 
choice of A5. For example for A5 = 0, the higgsino fraction fn (= l-ZV^sp + jA'^jeP) ^ denoting 
the index of the LSP, in this case 1) of the LSP decreases (see appendix [A]) thus suppressing the 

elastic scattering rate further by an order of magnitude. Of course = will increase the B — B' 
mass splitting, yet a splitting of only a few GeVs has little impact on the relic abundance(fig. [2)o). 

Conversely enhanced rates can be found for As > g' /\/2. 

The bino LSP up to a few hundred GeV is therefore a natural DM candidate in the DG model. 
The discovery of such neutralinos and sleptons of a few hundred GeV is within the reach of the 
LHC. Furthermore, only a rough estimate of the slepton/LSP mass difference would be sufficient to 
point out an inconsistency between the relic abundance of dark matter obtained from cosmological 
mesurements and the one predicted from collider measurements of the SUSY spectrum if done in 
the context of the MSSM. 
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6.1.3 Mixed bino/higgsino LSP 



A mixed bino/higgsino LSP is a more natural DM candidate than the bino because of the efficient 
annihilation into gauge boson pairs or top quark pairs. In particular the annihilation into gauge 
boson pairs proceeds through t-channel chargino exchange or Z/H exchange. For all annihilation 

diagrams some higgsino or W, W fraction of the LSP is involved, see the explicit expressions for 
the LSP couplings in Appendix C. Here and in the following we assume heavy sleptons Mr = 

Mj^ = 1 TeV as they are not needed for efficient annihilation. 

As a sample scenario we take M2D = 2Afi£) so that the wino fraction of the LSP is small and 
fix the mass of all sfermions to 1 TeV. In the DG model we find the contour of = 0.11 in 
the fi — Mid plane, see fig. El The contour corresponds to fi ^ Mir, and the region below the 
contour gives Q.h? < 0.11. The LSP is dominantly bino/bino' with a higgsino fraction that ranges 
from 2-30% along this contour as one increases the LSP mass. The small higgsino fraction of the 
LSP implies a small annihilation cross section into W pairs. This is however compensated by the 
contribution of the coannihilation channels to the relic abundance, with XiX^; 
into gauge bosons or heavy quarks. This is to be contrasted with the MSSM, where a similar 
contour in the — Mi plane is found but which feature a Higgsino fraction fn ~ 30%. Since 
coannihilation channels do not enter direct /indirect detection rates, the rates are suppressed in 
the DG model relative to the MSSM, compare fig. [3] with fig. [H The sharp variation in cr^Op at 

the edge of the figure correponds to the onset of the efficient annihilation of neutralinos through 
a Higgs resonance. Note that the mass spectra in the two models are rather similar, apart from 

the fact that there are two nearly degenerate B, B' and an additional chargino in the DG model. 
6.1.4 Mixed bino/wino/higgsino LSP 

We consider finally the case of a mixed bino/wino LSP, for example we choose M2D = l.lMiu- 
This choice for the gaugino masses necessarily implies a roughly 10% mass splitting between the 
LSP, an equal mixture of B/B', and the X2i^/^')^X3J Xii^, W') and the Xi ^xt ■ Thus naturally 
one finds important contributions from a variety of coannihilation channels. The contour of 
VLh? = 0.11 is displayed in fig. El Along this contour the coannihilation into gauge boson pairs or 
fermions dominate. Note that the main difference with the bino/higgsino case occurs for a light 
LSP, in particular, the small wino fraction is sufficient to provide efficient annihilation through 
chargino exchange, thus larger values of ^ are allowed. On the other hand in an MSSM model 
with a mixed wino LSP, here we choose M2 = l.lMi, only about half the effective annihilation 
cross section comes from coannihilation with Xi^xH^ dominant mode is LSP annihilation into 
W pairs (or top pairs). A smaller gaugino mass splitting, for example M2 = l.OSMi requires an 
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Figure 3: a) Contours of Jl/i^ = 0.11 in the n — Mid plane in the DG model for the mixed 
bino/Higgsino M2D = 2Mid (full blue line) and for M2D = l.lMiD(full green line), M2D = 
l.lMiD, Xs = l(red dashed line) b) Content of the LSP, fh (full lines) fw (dashed lines), same 
colour code function of m^o - same colour code as a) with the CDMS limit (full 

black line) d) cru |o as a function of m^o - same colour code as a) 
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Figure 4: a) Contours of = 0.11 in the // — Mi plane in the MSSM for the mixed bino/higgsino, 
M2 = 2Mi (blue) and for M2 = l.lMi (green) and M2 = l.OSMi (red) b) Content of the LSP, fh 
(full lines) fw (dashed lines) - same colour code as a) c) a^p as a function of m^o - same colour 
code as a) with the CDMS limit (full black line) d) crv\Q as a function of m^o - same colour code 
as a) 
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even smaller higgsino fraction. 

The direct detection rate which relies mainly on the LSP higgsino fraction is lower than for the 
mixed Higgsino LSP, this is especially true for light LSPs, see fig. [3^. As in other scenarios, the 
elastic scattering rate in the DG model is suppressed as compared to the MSSM. This statement 
is strongly dependent on the value of Xs- For example for = 1 we found a rate that increases 
by one or two orders of magnitude, this means within the range of the next run of Xenon |56j . 
The indirect rate is also increased in this case. 

6.2 Large Major ana masses Mi = M2 = 1 TeV 

First consider the M[ = M2 = case. When /i = 1 TeV, the LSP is dominantly B' with some B 
admixture and with a mass niLSP ^ ^Id/Mi. Both the wino and higgsino fraction of the LSP are 
very small, therefore the usually efficient annihilation channel into W pairs is suppressed. A value 
of O/i^ = 0.11 can only be reached because of the coannihilation channels such as X2X1') xtxi i^^to 
gauge boson pairs. These channels proceed through the wino component of the heavier neutralinos 
and charginos. For coannihilation to work one needs roughly a 10% mass difference between the 
LSP and the neutral and charged NLSP, NNLSP. In that sense the weak scale chargino/neutralino 
sector is similar to the one of the bino/wino scenarios in the MSSM, since additional states are 

at the TeV scale. In fig. [5] we display the contour Qh'^ = 0.11 in the m2D-fniD plane for different 
values of /x. When the higgsino fraction of the LSP increases, more precisely when the LSP mass 
becomes comparable to /i, annihilation becomes very efficient and the relic abundance is always 
nh'^ < 0.11. For u = 300 GeV, this occurs for Mid ~ 550 GeV, or m^o « 260 GeV. 

Al 

As discussed in previous scenarios, the predictions of the DG model for elastic scattering cross 
sections are usually suppressed when compared to an equivalent MSSM scenario. This is related 
to the fact that the LSP has a lower Higgsino fraction and that the relic abundance relies more 
heavily on coannihilation in the DG model. For example for /i = 1000 GeV the rate is suppressed 
by one order of magnitude in the DG model, see fig. [Sb- Here the MSSM rate corresponds to the 
Qh'^ = 0.11 contour for /x = 1000 GeV. As before, a large increase in the elastic scattering rate 
is expected when the LSP has a significant higgsino fraction. This occurs when Mlsp ~ or 

when Xs 7^ g' /\/2, see for example the contour /i = 1 TeV, As = 1 in fig. [Sb. Note that when the 
detection rate is small, interference between the squark and Higgs exchange can lead to a further 
suppression of the detection rate, see the dips for the /i = 300, 500 GeV scenarios in fig. [5)3. 

The self-annihilation of the LSP at v = is small for the B{B')/W scenario, the domimant 

channels are W^"^!^" or when this channel becomes kinematically accessible. As usual, the rate 
increases significantly with the higgsino content of the LSP, see fig. [5h. 

We will not discuss the more general case of the DG model where all Majorana and Dirac 
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Figure 5: a) Contours of Jl/i^ = 0.11 in the M2D — Mid plane in the DG model for /j, = 
1000 GeV(blue), /j = 500 GeV(green) and n = 300 GeV(red) b) cj^p as a function of m^o - 

same colour code as a). The dash curves correspond to fi = 1000 GeV (MSSM) and A5 = 1 (DG 
model) c) |o in the DG model as a function of m^o - same colour code as b). 

mass terms are present. When these masses are at the electroweak scale, there is no difficulty 
in finding dark matter scenarios with il/i^ ~ 0.1 and a LSP below the TeV scale. The LSP can 
be either a B or mixed B/h, B' /h, B /W{W'). All neutralinos and charginos could be within the 
kinematical reach of the LHC, the discovery of additional neutralino/chargino states would be the 
most obvious way to distinguish this model from the MSSM. 

7 Conclusion 

Dirac masses for gauginos can be obtained by pairing the MSSM gauginos with additional singlet, 
triplet and octet states in the adjoint representation. The model then contains additional neu- 
tralinos and charginos as well as new scalar particles. We have discussed both the case where the 
adjoint scalars decouple leaving a new quartic Higgs coupling in the effective potential and the 
one where the singlet remains light. This model has R-parity so the LSP is stable, the LSP could 
be a gravitino or a neutralino. 

We have made a first exploration of the parameter space of the model in the case of the neu- 
tralino LSP to find regions where the dark matter relic abundance is in agreement with cosmolog- 
ical measurements (assuming the standard cosmological scenario). Among the possible scenarios 
the ones that have a feature that distinguish them from the MSSM include the Dirac gaugino 
LSP that has a non suppressed annihilation into light fermions, the pseudo-Dirac bino LSP that 
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(co) annihilates into leptons via slepton exchange with sleptons much heavier than expected in 
the MSSM, as well as several mixed bino/wino/higgsino or bino/higgsino scenarios. For the latter 
scenarios the direct /indirect detection rates are often expected to be lower than in the MSSM, 
yet could be within the range of future detectors. In our numerical analysis we have concentrated 
on the DM candidates in the 100-500 GeV range, thus with several states within the range of 
the LHC. We have avoided a detailed discussion of scenarios where the annihilation of the LSP 
is made efficient by the presence of a Higgs resonance, as these scenarios require fine-tuning the 
masses of the LSP and Higgses. We expect that additional annihilation channels involving the 
new scalars in the model with a weak scale extra singlet would give new regions of parameter 
space with efficient DM annihilation as is found in the NMSSM. Finally we mention that the mass 
splitting between the two lightest neutralino states can be small although unless Xs is fixed exacty 
to its N=2 value the splitting is expected to be larger than the lOOkeV needed for inelastic dark 
matter scattering. 
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A Rotation Matrices and Eigenvalues for Bino/Wino Neutralinos 

In this appendix we give the first order rotation matrices for the neutralinos in the approximation 
that n is much larger than the other masses. Recall that the Lagrangian contains a term 



£ D --x^M^,XJ = -^X^A^S^'X^- (A.l) 



we write 



X'i = SiiRijXj = NijXj (A.2) 

where Rij is a real orthogonal matrix and Su a unitary diagonal matrix of phases that ensures that 
all of the masses are positive. For the case of that the neutralino is mostly bino, we assume that 
miD,m2D, Ml, M[, M2, M2 are of the same order and much smaller than /x. Then the eigenstates 
prior to mixing have mass eigenvalues 



± 1 
mi =2 



(Ml + M() ± J {Ml - M{)2 + Aml^ 



(A.3) 
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and similarly for mf. We then define 



ci = cos ffi = 



si = smf^i = 



miD 



so that 



-si ci M mii? Ml J \ si ci 



mj" 
m7 



(A.4) 



(A.5) 



and again similarly for the Wino states. Then Su = mf /\nii\, S22 = I and similarly 



for the mf values. The matrix Rij is given by 



i?ii=ci - 

Rl2 = Sl + 



Mlsisl,{2ciCf3SiSf}i{g'y + 2A|) + V2g'{cl - gf)(c| - 4)Ag) 
Mlcis^^{2cicpsisp{{g'f + 2A|) + V^g'icj - si){c} - s^^)Xs) 



R 



13- 



cwswMl 



{g'Y{m^ -mf)iJ. 



R 



gg'{mf — m2 ){'mf — mg") 
-ginq - m^)c2S2{2cpg'siSp + a/2cic|As - a/2cis|A5 



14 = 



(m^ - clrr4 - slm^ ){V2cl{g')si - V2{g')siSp - AciCf)SpXs)XT 



Rl5' 



gg'{mf - m^)(m]'" - mj) . 

-XT{m2 - m^)c2S2{V2cp{g')si - V2{g')siSp - 4C1C/3S/3A5)) 
Mzsw 



Mzsw 

-n-16 = ; 

91^ 



g'sisp + V2ciCfjXs 
g'sicp - V2cispXs 



(A.6) 
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^21= -Si - 
-R22 = Cl - 



R 



23- 



(m^ - m^){V2ciCp{g')XT - V2ci{g')sp\T + ^cpsispXsXr) 



R24-- 



R25 = 



g{g'){m^ - m2){m^ -rrq)li 
+(m^ - m^){-2ciC2Cpg{g')s2Sp + Ac^c^sispXsXr + V2{cp - sp){c2gsiS2Xs - cicl{g')XT) 

■ _ ^wSwMz^ _ ^+^2cici3g{g')sf3 + V2gsiXs - 2V2clgsiXs) 

g{g'){m^ -m^jim^ -mjjul 

-{rrq - m^)(V2(c| - s^f)){-gsislXs + ciC2{g')s2XT) + 2cpsp{cig{g')sl + 2c2SiS2XsXt) 

Mzsw " 



R26 = — 



91^ 

Mzsw 



R. 



31- 



g'cisp - -^sicpXs 
g'cicp + ■\/2siSi3Xs 



(A.7) 



{c\m^ + s^m| - m^)Xs{^/2c^pgs2 - V2gs2Sp - Ac2CpSpXT) 



R. 



32- 



R33 = C2 — 



g{g'){m'l - m^){m^ - )/i 
-^ci(g')(m^ - ml)si(2cpgs2Sp + \/2c2c|Ar - \/2c2s|At) 

+ + 37- ci(mi -m{)sxXs[y2cQgs2-\/2gs2Sa-^C2CjiSf)XT) 

g{g'){mj — mj){mj — mi)ij,i 

+{m^ - slm^ - clmf){-2cpg{g')s2Sf) - V2c2Cp{g')XT + V2c2ig')spXT) 

MZ^S2C^{2c2CpS2Sp{{g)^ + 2A2 ) + V2g{4 - sDic} - 4)At) 



R34 = S2 + 



(3)2 -rr4)ii 
MZ^C2^w{2c2CpS2S(}{{gf + 2A2,) + ^g{cl - sDic} - sI)Xt) 



(5)2(777,2 -nq)ii 



R35-- 



Mzcw 



R. 



36'- 



91^ 
9IJ- 



gs2S/3 + V2C2C13XT 

gs2Cp - V2C2S13XT 



(A.8) 
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R. 



■41 - 



ci(5')("^i - rn];)si{2c2Cf3gsij + \/2s2(-c| + 4)"^^) 



g{g'){mj - m2 )(m2 - m^)// 
+{cjm]; + slmj - rrq)Xs{V2c2g{cj3 - sp){cp + sp) + 4c/3S2S/3Ar) 
cwSwMl 



ci{m^ - mf)siXsiV2c2g{ci3 - Sf}){ci3 + s^) + 4c/3S2S/3At) 



i?43 = -S2 - 



g{g'){mj - )(m2 - )n 

+{g'){c\rnl + s^mj; - m2){2c2Cpgsi3 + a/2s2(-4 + s})Xt) 

MZ\2C^{2c2CfsS2Sp{{g)^ + 2A2.) + V2g{4 - - 4)At) 



Ru=C2 — 

Mzcw 



{gYim^ -mj)n 
MZ^S2cl^{2c2CfiS2Sp{{gf + 2X1) + ^5(ci - - 4)^^) 



(5)2 -mJ)iJ, 



i?45 = — 
i?46 = 



91^ 



gC2Sp - V2s2CpXT 
gC2Ci3 + ■\/2s2SpXT 



(A.9) 



-R51 

R53 

R54 

R55 
R56 



MzSwXsjcp - s/s) 

MzcwXricp - sp) 

Mzcwjcf) + Sfj) 
V2ti 

I/V2 
-I/V2 



(A.10) 
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-^61 
-^62 
-^63 



Mzsw>^s{ci3 + sp) 
Mzswjcp - sp) 

Mzcw>^t{ci3 + sp) 
Mzcw{cp - sp) 



V2m 



1/V2 



1/V2 



The second order eigenvalues for the above states are given by 



m 



+ 



c/3S;3((5')'s? - 2c?A|) + V2ci(5')si(4 - sl)\s 



= mi — 



c/3S/3((5')^sf - 2cf A|) - V2ci(5')si(4 - 4)A5 



2M2 s2 



4/^(5')' L 

-2 „2 



2M2f2 



,/^2 



5Ai 

(5')' + 2(c^-5^)2a| + 



2M2r2 



5^(c/3-s/3)^ + 2A^ 



(A.ll) 



(A.12) 
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B Rotation Matrices and Eigenvalues for Mostly Higgsino Neu- 
tralinos 



In this section we give the rotation matrix to first order for the case that /x is much less than the 
other masses. 

Rii = ci 

Ru = si 

Ri3 = 

Ru = 

Rib = ^^^T~^i'^^SCiSf}- gcpsi) 

Rm = ^^{V2Xscpci+g'spsi) (B.l) 



R21 = -si 

R22 = ci 

R23 = 

-R24 = 



R25 = - ^^^"^ {g'cpci + V2\sspsi) 

R26 = ^^^^^{g'cisis - V2XsC)3Si) (B.2) 
g'm^ 



R31 = 

R32 = 

R33 = C2 

Ru = S2 

R35 = ^^"y {-V2XTC2S3 + gcaS2) 

R36 = -—^^(A/2ArC/3C2 + 35/352) (B.3) 
gmj 
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R41 


= 


R42 


= 


R43 


= -S2 


Ru 


= C2 


R45 


Mzcw 







i?46 = — 



Mzcw 



{90/302 + a/2AtS/3S2) 

(-902813 + V2XTO13S2) 



(B.4) 



Mzsw 

2g' 771^171^ 

Mzsw 
2g'm^mf 

Mzcw 

2g' 1712 1712 

Mzcw 

2g' 7172 7172 



2\s{m^cl + 7n'lsl){ci3 - sp) + \/2{g'){7n^ - m^")ci(c^ + Sf3)si 
^{9){oi3 + si3){rnlcl + m^Si) + {m^ - 7n'^)Xs{op - sp)2cisi 
2\t{ci3 — S/3)(c|m2 + s^Tn^) + V2g{m2 — mJ)(c/3 + sp)c2S2 
V2g{slm2 + c^7n^){cp + sp) + 2AT(m2 — mJ)(c/3 — 5^3)0252 



1 



V2 



+ Rh 



(B.5) 



2g'7i7]^7nf 

Mzsw 
2g'7ni7nf 

Mzow 
2g'7n27n2 

Mzow ' 

2g' 7172 7112 



- 2{7nicl + 7n'lsl)Xs{of3 + sp) + V2{g'){7ni - 7n'l)ci{cf3 - sp)si 

- 2(m^ - 7i7f)Xsoi{c/3 + Si3)si + ■\/2{g){cp - S/3)(m]'"ci + m^^s?) 
2(cim2 + simJ)AT(c/3 + - V2g{7i72 - m^)(c;3 - 5^3)20252) 

- 2V2g{sl7n2 + C2m^)(c/3 - sp) + 2(m2' - m^)Ar(c/3 + Sf3)2c2S2 



(B.6) 
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where 
Rh 



8-\/2g'^{g'ym^ mfm2 

C2/3|(y)^'rtj"mf (-(m2 + m^){g'^ - 2X^) + (7713 - in^){g'^ + 2A|.)(c| - sl))c^ 
+g^m2mt{'{m^ + m+){{g'f - 2A|) + {m^ - mfWf + 2A|)(c? - s?))^^ 
+S2i39{9')'2V2i^{g')mimf{m2 - rr4)\TCy^2c2S2 + g{m{ - m^)m2 m^A52ciSis^ 
The non-LSP eigenvalues are 



(B.7) 



m^2 



m^3 



mf + 



mi + 



1712 + 



mo + 





2A|cf + (p')' 


^1 








2M2 s2 


cl{gr + 2slXl 
















9IJ' 




2M2r2 









(B.8) 
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while the LSP and NLSP are given by 



= /i 



2 

2 



2 

2 



M|g^(V2(g^)ci(c^ + sp) + 2A5(-c^ + sp)si) 
A{g'Ym^ 

Mlsl,{2\sci{cp - sp) + ^/2{g'){cp + 
4(5')2m^ 

Mlcl^{^/2gc2{cp - sp) + 2Xt{cp + sp)s2f 



Ag^m. 



M|c^(-2AtC2(c^ + sp) + V2g{cp - 8/3)32) 



2 

2 



2 

2 



4((7')^mj~ 

M|s2^(2A5Ci(c;3 + sp) + \/2(5')(-C/3 + sp)si)'' 



(B.9) 



4(5')2m]^ 

C Couplings of the LSP 

The neutrahno-chargino-W interactions depend on the W, W', h components of the neutraUno, 

= Ya^ {Cfil - 75) + (1 + 75)) ^W- + h.c. (C.l) 

where 

Cf = ^ (2Ni^Ua2 + 2N,sUal + V2iV,5C/a3) 

= ^ (2iV,4K2 + 2iVi3Ki - V2NieVa^') (C.2) 

When only gaugino Dirac mass are presents and when miD,m2D « fi the dominant contribution 
comes from the NisUn term. 

The Majorana neutralino coupling to the Z is driven by the higgsino component, as in the 
MSSM, 

C = l'^il''l5C'i^iZ^ + h.c. (C.3) 
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2cw 



(C.4) 



The neutrahno couphngs to the hght Higgs also depend on the W,W',h components of the 
neutrahno 



where 



swcw 



[e {cwNii - swNi2) {Zhj^^Ni5 - Zh^^Nio) 



(C.5) 



+ V2SWCW {NiiXs - NisXr)) {Zh.^N,^ + Zh^^Nis) 



(C.6) 



Here Z^, is the scalar mixing matrix. 
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